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摘要 :【 目 的] 本 研究 旨 在 探讨 调控 件 裙 追 寄 蝇 Exorista civilis 滞 育 的 重要 关联 基因 以 及 代谢 途径 ， 
为 明确 该 虫 在 转录 组 水 平 下 的 滞 育 分 子 机 制 提供 理论 基础 。[【 方 法 ) 采 用 新 一 代 高 通 量 测 序 平台 
Ilumina HiSeq™2000 对 伞 裙 追 寄 晶 非 滞 育 以 及 滞 育 的 肾 进 行 转录 组 测序 分 析 以 及 生物 信息 学 分 
析 ; 应 用 KAAS 在 线 pathway 比 对 分 析 工 具 对 筛选 出 的 满足 padj <0. 05 且 fold change 宇 32 或 padj < 
0.05 且 fold change 友 1/32 条 件 的 差异 表达 基因 进行 KEGG 通路 富 集 分 析 。【 结 果 】 根 据 测 序 结果 ， 
共 获 取 58 050 个 基因 。 差 异 倍数 在 32 倍 以 上 的 滞 育 关联 基因 (diapause-associated genes, DAGs) 
有 454 个 ,其 中 406 个 表达 上 调 , 共 涉及 134 条 通路 ,包括 氧化 磷酸 化 、 柠 榜 酸 循环 及 其 他 重要 通 
路 ;48 个 表达 下 调 , 涉 及 32 条 通路 。KEGG 通路 富 集 分 析 结 果 表 明 , 滞 育 关联 基因 主要 集中 在 信号 
转 导 、 内 分 沁 系 统 、 碳 水 化 合 物 代谢 等 途径 中 。[【 结论 ] 本 研究 获得 的 伞 裙 追 寄 蝇 转录 组 数据 揭示 
了 其 滞 育 调控 的 重要 代谢 途径 与 关联 基因 。 

关键 词 : 伞 裙 追 寄 晶 ; 转录 组 测序 ; 滞 育 ; 滞 育 关联 基因 ; KEGG 分 析 
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Abstract: [Aim) This study aims to explore the crucial associated genes and metabolic pathways that 
regulate the diapause of Exorista civilis, and to provide a theoretical basis for clarifying the molecular 
mechanism of diapause of this species at the transcriptome level. 【Methods】The Ilumina HiSeq™ 2000 
platform was used to perform transcriptome sequencing and bioinformatics analysis of the non-diapause 
and diapause pupae of E. civilis. KAAS ( KEGG Automatic Annotation Server) online pathway 
comparison tool was used to analyze the KEGG pathway enrichment of the differentially expressed genes 
that meet the conditions of padj <0.05 and fold change 三 32 or padj <0.05 and fold change < 1/32. 
【Results】According to the sequencing results, a total of 58 050 unigenes were obtained. There were 454 
diapause-associated genes (DAGs) with the difference of more than 32-fold changes. Among them, 406 
DAGs were up-regulated and involved in a total of 134 pathways, including oxidative phosphorylation, 
citrate cycle and other important pathways, and 48 DAGs were down-regulated and related to 32 





























基金 项 目 : 国家 重点 研发 计划 (2017YFD0201000, 2017YFE0104900) ; 内 蒙古 自治 区 自然 科学 基金 项 目 (2019BS03011) 

作者 简介 : 张 博 , 女 , 1994 年 11 月 生 , 内 蒙古 乌 兰 察 布 人 , 硕士 研究 生 , 研究 方向 为 农业 昆虫 与 害虫 生物 防治 ，E-mail: zhangbojenny1001@ 
163. com 
* 通讯 作者 Corresponding author, E-mail: liuaiping806@ sohu. com 
收 稿 日 期 Received: 2019-12-08; 接受 日 期 Accepted: 2020-06-01 






























































8 期 张 ” 博 等 : 伞 裙 追 寄 晶 滞 育 关联 基因 的 转录 组 学 分 析 933 











pathways. KEGG enriched pathways revealed that the DAGs are mainly involved in signal transduction, 


endocrine system and carbohydrate metabolism. 【 Conclusion 】 The transcriptome data of E. civilis 


obtained in this study have revealed the important metabolic pathways and associated genes in diapause 


regulation in E. civilis. 


Key words : Exorista civilis; transcriptome sequencing; diapause; diapause-associated genes; KEGG analysis 





伴 裙 追 寄 蝇 Exorista civilis 是 一 种 牧草 优势 寄生 
性 天 敌 ,隶属 于 双 翅 目 ( Diptera) 寄 蝇 科 (Tachinidae ) ， 
能 够 对 粘 虫 朵 ythimna Separata、 草地 量 Loxostege 
sticticalis 、 草 地 贪 夜 峨 Spodopiera frugiperda 等 多 种 牧 
草 害虫 起 到 有 效 的 防 控 作用 (王建 梅 等 ,2013， 
2015 ) 。 昆 虫 滞 育 是 一 种 生理 过 程 , 其 特征 在 于 行 
为 修饰 .代谢 抑制 以 及 增强 的 胁迫 耐 受 性 (King and 
Macrae, 2015 ) 。 近 年 来 ,关于 伞 裙 追 寄 蝇 滞 育 类 
型 洲 育 温 光 周期 的 诱导 以 及 灌 育 关联 蛋白 的 人 研究 
取得 了 一 定 的 成 果 ( 韩 海底 等 , 2018; 姜 珊 , 2018 ) ， 
为 揭示 其 滞 育 的 分 子 机 制 黄 定 了 基础 。 

新 一 代 测 序 技术 的 出 现 , 极 大 促进 了 昆虫 转录 
组 学 研究 ( 张 棋 变 和 豆 明 龙 , 2013 ) 。 目 前 利用 高 通 
量 技术 对 双 埃 目 昆虫 的 滞 育 研究 也 取得 了 一 定 的 进 
展 , 利 用 高 通 量 技术 发 现 了 在 柑橘 大 实 晶 Bactrocera 
minax 清 育 调控 的 相关 研究 中 筛选 出 12 个 显著 上 
调和 6 个 显著 下 调 的 滞 育 关联 基因 以 及 重要 的 清 育 
关联 和 蛋白 hsp23 ，hsp70 和 hsp90 (Lu et al.,，2014; 
Wang et al.，2019 ) 。 张 倩 等 (2019 ) 利用 iTRAQ 技 
术 筛 选 淡色 库 蚊 Culex pipiens pallens 越冬 后 期 和 沸 
育 期 的 差异 表达 和 蛋白, 共 鉴定 出 244 个 差异 表达 和 蛋 
白 , 包 括 126 个 上 调和 蛋白 和 118 个 下 调和 蛋白 ,这 些 差 
异 表 达 和 蛋白 的 功能 与 脂 质 代谢 、 细 胞 骨架 重 塑 .碳水 
化 合 物 代谢 等 有 关 。 

目前 ,调控 伞 裙 追 寄 蝇 滞 育 的 分 子 调控 机 制 仍 
然 不 完全 明确 。 为 此 ,本 研究 对 伴 裙 追 寄 蝇 正 常 发 
育 及 滞 育 的 师 进 行 转录 组 测序 ,对 筛选 出 的 清 育 关 
联 基因 进行 KEGG 通路 富 集 分 析 , 了解 灌 育 关联 基 
因 主 要 参与 的 代谢 途径 和 信号 转 导 机 制 ,为 进一步 
研究 伞 裙 人 妃 寄 晶 的 滞 育 调控 基因 以 及 分 子 机 理 疯 定 
基础 ,并 为 通过 遗传 信息 研发 滞 育 剂 . 灌 育 解除 剂 ， 
从 而 实现 人 工 调控 锌 裙 追 寄 蝇 等 天 敌 昆 虫 的 寿命 ， 
实现 在 田间 精准 释放 ,发 挥 生 物 防治 的 最 佳 效 果 提 
供 科学 依据 。 


1 材料 与 方法 


1.1 供 试 昆虫 及 滞 育 诱导 
伞 裙 追 寄 蝇 贤 为 中 国 农业 科学 院 草 原 人 研究 所 实 







































































































































































验 室 所 饲养 。 将 正 党 发 育 的 非 汪 育 晴 在 25Y +1%C 
室温 条 件 下 培养 1 d, 经 用 液 所 处 理 后 放 和 人 -80Y 冰 
箱 保存 备用 ;将 正常 发 育 的 非 滞 育 贤 置 于 人 工 气 候 
箱 内 ,在 温度 11 + 1% 、 相 对 湿度 70% + 上 10% 、 光 周 
期 8L: 16D 的 条 件 下 ,经 过 40 d 诱导 可 获得 清 育 态 
晴 ( 姜 珊 , 2018) ,40 d 后 选取 重量 与 正常 发 育 非 滞 
育 晴 一 致 的 存活 济 育 师 , 经 用 液 氮 处 理 后 放 入 
-80% 冰箱 保存 备用 。 然 后 ,将 备用 的 伞 裙 追 寄 蝇 
正常 发 育 非 沛 育 肾 和 沛 育 肾 进行 RNA 样品 检测 , 主 
要 方法 为 :(1) 琼 脂 糖 凝 胶 电泳 分 析 RNA 降解 程度 
以 及 是 否 有 污染 ; (2) NanoDrop 检测 RNA 的 纯度 
(0OD;aZOD:m 比 值 );(3)Qubit 对 RNA 浓度 进行 精 
确定 量 ;(4) Agilent 2100 精确 检测 RNA 的 完整 性 。 
1.2 转录 组 测序 和 组 装 

伞 裙 追 寄 蝇 非 滞 育 和 淆 育 晴 的 转录 组 测序 工作 
由 北京 诺 禾 致 源 科 技 股份 有 限 公 司 完成 。 测 序 进 行 
3 个 生物 学 重复 ,每 个 生物 学 重复 以 4 头 晴 为 样本 ， 
两 种 晴 共 构建 6 个 文库 。 样 品 检测 合格 后 ,用 带 有 
oligo(dT) 的 磁 珠 富 集 真 核 生 物 mRNA。 随 后 加 入 
fragmentation buffer 将 mRNA 打 断 成 短片 段 , 以 mRNA 
为 模板 ,用 六 碱 基 随 机 引物 (random hexamers ) 合成 
cDNA 第 1 链 , 然 后 加 入 缓冲 液 ,再 加 纯化 的 双 链 
cDNA。 纯 化 的 双 链 cDNA 先进 行 末端 修复 加 和 A 尾 
并 连接 测序 接头 ,再 进行 片段 大 小 选择 。 最 后 进行 
PCR 扩 增 ,得 到 最 终 的 文库 。 文 库 构建 完成 后 , 先 
使 用 Qubi2. 0 进行 初步 定量 ,稀释 文库 至 1.5 ng/ 
phEL ,随后 使 用 Agilent 2100 对 文库 的 insert size 进行 
检测 ,insert size 符合 预期 后 ,使 用 Q-PCR 方法 对 文 
库 的 有 效 浓度 进行 准确 定量 (文库 有 效 浓度 >2 
nmolL) ,以 保证 文库 质量 。 库 检 合 格 后 ,把 不 同文 
库 按照 有 效 浓度 及 目标 下 机 数据 量 的 需求 pooling 
后 进行 高 通 量 测序 。 测 序 完成 后 , 高 通 量 测序 
( lumina HiSeq”2000) 得 到 的 原始 图 像 数 据 文件 
经 碱 基 识 别 (base calling) 分 析 转 化 为 原始 测序 序 
列 ,然后 进行 测序 错误 率 分 布 检 查 ,获得 高 质量 的 序 
列 (clean reads ) 后 ,采用 Trinity 软件 对 clean reads 
进行 拼接 ,Corset 参差 聚 类 后 得 到 最 长 的 序列 进行 
后 续 法 分 析 。 
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1.3 生物 信息 学 分 析 

为 获得 全 面 的 基因 功能 信息 ,对 伴 裙 追 寄 蝇 进 
行 了 七 大 数据 库 的 基因 功能 注释 ,包括 :NCBI 非 元 
余 蛋 白 序列 数据 库 (NCBI non-redundant protein 
sequences,，Nr ) , NCBI 非 宛 余 核 背 酸 序列 数据 库 
(NCBI non-redundant nucleotide sequences，Nt) , 蛋白 
家 族 ( protein family，Pfam ), 直系 同 源 基因 徐 
(cluster of orthologous groups of proteins, KOG/ 
COG ) ，Swiss-Prot 数据 库 (a manually annotated and 
reviewed protein sequence database ) ， 京 都 基因 与 基 
因 组 百科 全 书 (KEGG ortholog database，K0O), 基 
本 体 论 (gene ontolooy, GO) 。 拼 接 得 到 的 爹 裙 追 寄 
晶 转 录 组 测序 结果 与 UniProt 或 者 Nr 数据 库 进行 
Blastx 比 对 (将 核 昔 酸 序列 翻译 为 蛋白 ,再 进行 比 
对 ) ,筛选 条 件 E-value < le ,得 到 注释 基因 。 通 过 
RSEM 软件 得 到 每 个 样品 比 对 到 每 个 基因 上 的 
readcount 数目 ,并 对 其 进行 FPKM (expected number 
of fragments per kilobase of transcript sequence per 
millions base pairs sequenced ) 转换 ,分 析 基 因 的 表达 
水 平 ;对 转录 组 测序 的 非 滞 育 师 和 清 育 师 的 样品 进 
行 了 样品 间 的 相关 性 检查 ,采用 R 语言 进行 Pearson 
相关 系数 的 计算 ;将 正常 非 滞 育 肾 和 浪 育 肾 的 转录 
组 测序 结果 进行 比较 ,采用 校正 后 的 P 值 (Pp- 
adjusted ，padj ) 与 差异 表达 倍数 (fold change, FC) 
方法 进行 差异 表达 基因 挑选 ,padj 是 对 PP 值 进行 了 



















































































多 重 假设 检验 ,能 更 好 地 控制 假 阳 性 率 。 对 筛选 出 
的 差异 表达 基因 进行 聚 类 分 析 , 使 用 聚 类 软件 包 
pheatmap ,针对 差异 表达 基因 的 并 集 , 以 基因 的 相对 
表达 水 平 值 进行 聚 类 。 和 采用 相应 的 距离 算法 计算 每 
个 基因 之 间 的 距离 ,然后 通过 反复 迭代 ,计算 基因 之 
间 的 相对 距离 ,最 后 根据 基因 的 相对 距离 远近 分 成 
不 同 的 subeluster, 从 而 实现 聚 类 。 对 伴 裙 追 寄 蝇 非洲 
育 晴 和 淆 育 晴 差 异 表 达 基 因 的 筛选 条 件 为 : padj < 
0. 05 且 fold change 三 32 或 padqj <0. 05 且 fold change= 
1/32, 与 非 渍 育 肾 的 相 比 , 灌 育 肾 基 因 差 异 表达 倍数 
在 32 倍 以 上 的 淖 育 关联 基因 ,并 进行 重点 分 析 。 
1.4 滞 育 关联 基因 KEGG 通路 分 析 

将 灌 育 关联 基因 序列 制 成 FASTA 格式 文件 ,应 
用 KAAS( KEGG Automatic Annotation Server. http: 
/www. genome. jp/kegg/ ) 在 线 pathway 比 对 分 析 工 
具 进 行 KEGG 通路 富 集 分 析 。 







































































2 结果 


2.1 伞 裙 追 寄 蝇 非 灌 育 和 滞 育 蚂 转 录 组 测序 结果 

本 研究 对 伞 裙 追 寄 蝇 非 滞 育 和 滞 育 是 进 行 转录 
组 测序 ,高 质量 样品 数据 超过 7.42 Cb ,每 个 样品 的 
Q20 均 大 于 96.08% ,Q30 均 大 于 90.01% ( 表 1) , 通 
过 转录 组 数据 结果 发 现 其 准确 度 较 高 , 可 用 于 后 续 
分 析 。 






































表 1 伞 裙 追 寄 蝇 非 滞 育 和 灌 育 肾 转 录 组 测序 数据 质量 评估 统计 


Table 1 Quality statistics of sequencing data of the transcriptomes of non-diapause and diapause pupae of Exorista civilis 





样品 原始 reads 数 高 质量 reads 数 数据 大 小 (Gb) 碱 基 错 误 率 ( % ) Q20 Q30 (G+C) 
Samples Numberof raw reads Number of clean reads Size of clean bases Error rate of bases (%) (%) % 
D1 61251702 59 781 232 8.97 0.03 97.18 92.26 37.21 
D2 51135130 49 489 976 7.42 0.03 96. 08 90.01 36.56 
D3 57222216 55 774 338 8.37 0.03 97.07 92. 06 37.07 
ND1 51607212 50 396 094 7.56 0.03 97.22 92.41 38.59 
ND2 54932926 53 057 940 7.96 0.03 97.03 92.01 38.78 
ND3 69990434 68 331 496 10. 25 0.03 96. 87 91.66 32:77 


示 质 量 值 20 和 30 的 碱 基 所 占 百分比 Percentages of bases with the quality value 20 and 30 ，respectively- 

样品 间 基 因 表达 水 平 相关 性 是 检验 实验 可 靠 性 ”表明 灌 育 晴 样 品 与 非 滞 育 晴 样 品 基 因 表 达 模 式 存 在 
和 样本 选择 是 否 合理 的 重要 指标 。 样 品 间 基因 表达 ” 明显 差异 。 
水 平 的 相关 性 反映 了 样品 间 的 相似 程度 , 即 不 同 处 差异 表达 基因 聚 类 分 析 用 于 判断 不 同 实验 条 件 
理 或 组 织 的 样品 在 基因 表达 水 平方 面 的 相似 度 。 相 。 ”下 差异 基因 表达 量 的 聚 类 模式 ( 杨 春 梅 等 , 2003 ) 。 
关系 数 越 接近 1 ,表明 样品 之 间 基 因 表 达 模 式 的 相 。 图 2 红色 表示 高 表达 , 蓝 色 表示 低 表 达 。 本 研究 检 
似 度 越 高 ,样品 间 的 差异 表达 基因 也 越 少 ( 张 宇 镭 ，” 测 结果 中 蓝 色 和 红色 部 分 能 够 明显 区 分 ,说 明 样 品 
等 , 2005) 。 本 研究 中 生物 学 重复 间 的 样品 的 相关 ”整体 的 测序 质量 好 。 
系数 大 于 生物 学 重复 外 的 样品 的 相关 系数 (图 1)， 
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图 1 人 镍 裙 追 寄 晶 非 滞 育 和 滞 育 师 转 录 组 样品 间 基 因 表 达 水 平 Pearson 相关 系数 热 图 


Fig. 1 Heatmap of Pearson correlation of gene expression levels between transcriptome samples of non-diapause 





and diapause pupae of Exorista civilis 
R? . Pearson 相关 系数 的 平方 Square of Pearson correlation coefficient. 





D2 DI D D3 ND3 ND2 ND NDIl 


样品 Samples 


图 2 候 裙 追 寄 蝇 非 滞 育 和 滞 育 肾 转 录 组 中 差异 表达 基因 聚 类 图 
Fig. 2 Cluster analysis of differentially expressed genes in the transcriptome of non-diapause and diapause pupae of Exorista civilis 
D: 滞 育 肾 混 合 样 品 Mixed samples of diapause pupae; ND: 非 滞 育 晴 混 合 样品 Mixed samples of non-diapause pupae. Z 分 数 是 一 种 对 数据 进行 标 
准 化 的 方法 ,其 计算 方法 为 减 去 均值 除 以 标准 差 ,均一 化 处 理 后 的 数值 在 -2 ~2 之 间 。Z-score is a method of standardizing data. Its calculation 


method is to subtract the mean value and divide by the standard deviation. The value after normalization is between -2 and 2. 
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2.2 滞 育 关联 基因 KEGG 通路 分 析 途径 。 与 企 裙 追 寄 晶 非 沾 育 晴 相 比 , 滞 育 晴 转录 组 


伞 裙 追 寄 蝇 非 滞 育 师 和 滞 育 师 转 录 组 测序 在 
KEGG 直系 同 源 系 统 (KEGG orthdogy，KO ) 中 成 功 
注释 了 6 643 个 基因 ,KO 系统 将 各 个 KEGG 注释 系 
统 联系 在 一 起 ,可 完成 转录 组 的 功能 注释 。 根 据 篇 
选 条 件 :padj <0.05 日 fold change 二 2 或 padj <0. 05 
且 fold change 志 0.5 ,获得 在 伞 裙 追 寄 晶 非 滞 育 师 和 
潍 育 肾 转 录 组 间 差 异 2 倍 以 上 的 差异 表达 基因 共有 
2 648 个 , 清 育 师 中 上 调 表 达 基 因 1 518 个 ,下 调 表 
达 基 因 1 130 个 。 对 这 些 差异 表达 基因 进行 KECC 
分 类 (图 3) ,所 属 KEGG 代谢 通路 分 为 5 类 ,包括 细 
胞 过 程 、 环 境 信 息 处 理 、 遗 传 信息 处 理 、 代 谢 
(metabolism) 和 有 机 系统 ,这 些 基因 主要 集中 在 信 
号 转 导 、 内 分 泌 系统 、 细 胞 群落 和 碳水 化 合 物 代谢 等 


细胞 生长 与 死亡 Cell growth and death 

细胞 运动 Cell motility 

细胞 群落 Cellular community 

运输 和 分 解 代谢 Transport and catabolism 

膜 运输 Membrane transport 

祝 号 分 子 与 相互 作用 Signaling molecules and interaction 
信号 转 导 Signal transduction 

折 和 去 、 分 类 和 降解 Folding, sorting and degradation 

复制 和 修复 Replication and repair 

转录 Transcription 

翻译 Translation 

聚 糖 的 生物 合成 与 代谢 Glycan biosynthesis and metabolism 
核 音 酸 代谢 Nucleotide metabolism 

其 他 氨基 酸 的 代谢 Metabolism of other amino acids 

辅助 因子 和 维生素 的 代谢 Metabolism of cofactors and vitamins 
莅 类 和 聚 酮 化 合 物 的 代谢 Metabolism of terpenoids and polyketides 
脂 代谢 Lipid metabolism 

纵览 Overview 

能 量 代 谢 Energy metabolism 

碳水 化 合 物 代谢 Carbohydrate metabolism 
氨基 酸 代谢 Amino acid metabolism 

异 生物 素 的 生物 降解 与 代谢 Xenobiotics biodegradation and metabolism 
发 育 Development 

消化 系统 Digestive system 

排泄 系统 Excretory System 

感觉 系统 Sensory system 

循环 系统 Circulatory system 

神经 系统 Nervous system 

免疫 系统 Immune system 

内 分 泌 系 统 Endocrine system 

环境 适应 Environmental adaptation 


中 差异 在 32 倍 以 上 的 浪 育 关联 基因 为 454 个 ,可 分 
为 代谢 .环境 信息 处 理 、 遗 传 信息 处 理 、 细 胞 过 程 和 
有 机 系统 五 大 类 ,每 类 中 滞 育 关联 基因 参与 的 通路 
占 比 前 3 的 如 表 2 所 示 。 其 中 上 调 基 因为 406 个 ， 
下 调 基 因为 48 个 。 通 过 KEGG 数据 库 分 析 , 发 现 
454 个 滞 育 关联 基因 共 参 与 到 141 个 通路 ,其 中 参 
与 焦点 粘连 通路 的 清 育 关联 基因 最 多 ,为 13 个 ;其 
次 是 参与 碳 代谢 通路 的 ,为 12 个。 

2.2.1 代谢 :根据 KO 分 类 结果 , 涡 育 关联 基因 共 
参与 代谢 分 类 中 的 碳水 化 合 物 代谢 ( carbohydrate 
metabolism ) ,能量 代谢 (energy metabolism ) 、 脂 代谢 
( lipid metabolism )、 核 苷 酸 代 谢 ( nucleotide 
metabolism ) 、 氨 基 酸 代谢 (amino acid metabolism ) 等 


EEE 12 
N° 
ES 33 





20 30 40 50 60 70 


Unigenes 数 量 Number of unigenes 


图 3 人 镍 裙 追 寄 晶 滞 育 师 转 录 组 中 与 非 滞 育 晴 相 比 差 异 表 达 基 因 的 KEGG 功能 分 类 


Fig. 3 KEGG function classification of differentially expressed genes in the transcriptome of 


the diapause pupae compared with non-diapause pupae of Exorista civilis 
将 基因 根据 参与 的 KEGG 代谢 通路 分 为 5 个 分 支 : 细 胞 过 程 (A) ,环境 信息 处 理 (B) ,遗传 信息 处 理 (C) ,代谢 (D) 和 有 机 系统 ( 卫 ) 。The genes 


were divided into five branches according to the involved KEGG metabolic pathway: cellular process (A ) environmental information processing (B) ， 








genetic information processing (C) metabolism (D ) and organic system (E). 
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表 2 企 裙 追 寄 蝇 滞 育 肾 转 录 组 中 滞 育 关联 基因 KEGG 通路 分 析 


Table2 KEGG pathway analysis of diapause-associated genes in the transcriptome of diapause pupae of Exorista civilis 































































































代谢 通路 滞 育 关联 基因 参与 数量 占 比 (%) 通路 ID 
Metabolic pathways Number of diapause-associated genes Proportion Pathway ID 
代谢 Metabolism 

碳 代谢 Carbon metabolism 12 2.6 ko01200 

氧化 磷酸 化 Oxidative phosphorylation 9 1.9 ko00190 

柠檬 酸 循环 Citrate cycle (TCA cycle) 8 1.7 ko00020 
有 机 系统 Organic system 

胃酸 分 泌 Gastric acid secretion 10 :2 ko04971 

卵 母 细胞 减 数 分 裂 Oocyte meiosis 10 2 ko04114 

小 板 活 化 Platelet activation 10 2 ko04611 
状 腺 激素 信和 号 通路 Thyroid hormone signaling pathway 10 2.2 ko04919 

细胞 过 程 Cellullar process 

焦点 粘连 Focal adhesion 13 2.8 ko04510 

紧密 连接 Tight junction 9 1.9 ko04530 

肌 动 蛋白 细胞 骨架 的 调控 Regulation of actin cytoskeleton 9 1.9 ko04810 
环境 信息 处 理 Environmental information processing 

Hippo 信和 号 通路 Hippo signaling pathway 8 1.7 ko04390 

Rapl 信号 途径 Rapl signaling pathway 7 1.5 ko04015 

Hippo 信和 号 通路 - 蝇 Hippo signaling pathway-fly 6 1.3 ko04391 
遗传 信息 处 理 Genetic information processing 

内 质 网 蛋白 质 加 工 Protein processing in endoplasmic reticulum 2 0.4 ko04141 

泛 素 介 导 的 蛋白 质 水 解 Ubiquitin-mediated proteolysis 2 0.4 ko04120 





11 个 途径 ,涉及 47 个 通路 。 

分 析 结 果 显 示 , 在 糖 酵 解 / 糖 异 生 ( glycolysis/ 
gluconeogenesis) 通路 中 , 丙酮 酸 激 酶 ( pyruvate 
kinase, PK ) 基因 表达 上 调 101. 60 倍 ,在 柠檬 酸 循环 
(citrate cycle, TCA cycle) 中 ,丙酮 酸 脱氧 酶 的 El 组 
分 (pyruvate dehydrogenase El component, PDHA) a 
亚 基 基因 表达 上 调 151. 49 倍 ,PK 是 糖 醇 解 过 程 中 
发 生 第 二 次 底 物 磷酸 化 的 关键 酶 ,催化 反应 生成 丙 
酮 酸 。 在 有 氧 条 件 下 , 丙酮 酸 脱氧 酶 复合 体 
(pyruvate dehydrogenase complex，PDHC ) 催化 由 糖 
酵 解 产生 的 丙酮 酸 进行 氧化 脱羧 反应 ,促进 乙酰 畏 
酶 A 的 生成 ,进入 柠檬 酸 循环 。 在 脂 代谢 中 , 脂 酰 
辅酶 A 还 原 酶 (fatty acyl-CoA reductase, FAR ) 基 
表达 上 调 , 上 调 倍数 最 高 可 达 146. 49 倍 。 在 能 量 代 
谢 中 ,参与 氧化 磷酸 化 的 滞 育 关联 基因 有 9 个 ,氧化 
磷酸 化 是 一 个 电子 传递 过 程 , 对 线粒体 呼吸 电子 链 
的 传递 十 分 重要 。 细 胞 色素 C 氧化 酶 是 线粒体 中 
电子 传递 链 的 末端 酶 ,催化 电子 从 细胞 色素 C 转移 
到 氧 。 细 胞 色素 氧化 酶 的 5 个 亚 基 基因 表达 上 调 ， 
其 中 细胞 色素 C 氧化 酶 亚 基 4 (cytochrome c oxidase 
subunit 4 ，Cox4 ) 是 1 个 关键 亚 基 。 在 昆虫 激素 的 生 
物 合 成 (insect hormone biosynthesis ) 通路 中 , 保 幼 激 














































































































素 酯 酶 (juvenile-hormone esterase，JHE ) 基因 表达 下 
调 , 下 调 了 82.28 倍 。 
2.2.2 有 机 系统 :根据 KO 注释 结果 ,有 机 系统 中 
的 滞 育 关联 基因 参与 免疫 系统 .内分泌 系统 、 循 环 系 
统 .排泄 系统 ,神经 系统 .感觉 系统 以 及 消化 系统 等 
7 个 系统 ,其 中 有 58 个 清 育 关联 基因 参与 内 分 泌 系 
统 ,涉及 内 分 泌 系 统 的 15 个 通路 。 

钙 调 蛋白 (calmodulin，CALM) 基因 上 调 表达 主 
要 涉及 内 分 泌 系 统 , 钙 调 蛋白 本 身 无 任何 酶 活性 ,可 
以 通过 结合 Ca ' 形成 Ca * -CaM 复合 体 ,调控 细胞 
代谢 过 程 以 及 下 游 的 靶 和 蛋白 、 靶 酶 。Hsp70 的 一 种 : 
heat shock 70kDa protein 1/8( HSPA1s) 在 内 分 泌 系 
统 的 雌 激 素 信 号 通路 (estrogen signaling pathway ) 、 
免疫 系统 的 抗原 处 理 和 展示 (antigen processing and 
presentation) 通路 中 既 有 上 调 也 有 下 调 , 上 调 倍数 最 
高 可 达 98.53 倍 ,下 调 倍 数 为 2.56 倍 。 
2.2.3 细胞 过 程 :根据 KO 注释 结果 ,参与 细胞 过 
程 的 73 个 济 育 关联 基因 主要 涉及 焦点 粘连 (focal 
adhesion) 、 肌 动 蛋 白细胞 骨架 的 调节 (regulation of 
actin cytoskeleton ) 、 紧 密 连 接 (tight junction ) 、 颖 际 连 
接 (adherens junction ) 、 吞 哈 体 (phagosome) 、 溶 酶 体 
(lysosome) \ 细 胞 凋 亡 ( apoptosis) 、 卵 母 细胞 减 数 分 
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裂 (oocyte meiosis) 间隙 连接 (gap junction ) 、 过 氧化 
物 酶 体 ( peroxisome ) 、p53 信号 通路 (p53 signaling 
pathway)、 细 胞 周期 (cell cycle)、 胞 吞 作 用 
(endocytosis ) 这 13 个 通路 。 

在 吞噬 体 通路 中 肌 动 蛋白 基因 表达 占 比 最 高 ， 
共有 9 个 基因 被 注释 为 肌 动 蛋白 B/yl (actin beta/ 
gamma 1，ACTB_G1) , 肌 动 蛋白 是 真 核 生物 中 最 丰 
富 的 蛋白 质 之 一 ,对 维持 细胞 形态 具有 重要 的 作用 。 
在 溶 酶 体 中 组 织 蛋 白 酶 表达 占 比 最 高 ,共有 6 个 基 
因 分 别 被 注释 为 组 织 蛋 白 酶 G (cathepsin C， 
CTSG) .组 织 蛋 日 酶 S(cathepsin S，CTSS) 、 组 织 
白 酶 K(cathepsin K, CTSK ) 。 

在 细胞 凋 亡 通路 中 ,细胞 色素 C( cytochrome ec， 
CYC) 基 因 表 达 上 调 ,CYC 介 导 呼吸 链 中 的 电子 传 
递 ,影响 ATP 的 产生 ,其 释放 可 诱导 细胞 凋 亡 。 在 
过 氧化 氨 物 酶 体 通路 中 ,过 氧化 氢 酶 (catalase， 
CAT) 基因 表达 上 调 , 上 调 123. 17 倍 。 

2.2.4 环境 信息 处 理 : 根 据 滞 育 关联 基因 的 KO 比 
对 结果 , 环境 信息 处 理 共 涉 及 到 信号 转 导 (signal 
transduction)、 信 号 分 子 及 其 相互 作用 (signaling 
molecules and interaction ) 、 膜 转 运 ( membrane 
transport)3 个 途径 。60 个 清 育 关联 基因 共 参 与 环 磷 
酸 腺 苷 (cAMP signaling pathway ) , rapl 信号 通路 
(rapl signaling pathway ) , MAPK 信号 通路 ( MAPK 
signaling pathway ) , 钙 信 号 通路 (calcium signaling 
pathway ) ，cGMP-PKG 信号 通路 ( cGMP-PKG 
signaling pathway ) , hippo 信号 通路 (hippo signaling 
pathway ) 等 在 内 的 19 个 通路 。 在 ECM- 受 体 相 互 作 
用 通路 中 , 肌 腿 蛋白 (tenascin, TN) 基因 表达 下 调 ， 
下 调 倍数 最 高 可 达 51. 98 倍 。 
2.2.5 遗传 信息 处 理 : 根 据 KO 注释 结果 ,在 遗传 
言 息 处 理 第 二 层次 的 分 类 中 涉及 复制 与 修复 
(replication and repair) , 折 著 分 类 和 降解 (folding， 
sorting and degradation ) , 翻译 (translation ) 和 转录 
(transcription )4 个 途径 ,第 3 层次 分 类 中 涉及 内 质 
网 蛋白 加 工 (protein processing in endoplasmic 
reticulum) 、RNA 运输 (RNA transport)、 泛 素 介 导 的 
蛋白质 水 解 (ubiquitin-mediated proteolysis ) 等 7 个 通 
路 。 泛 素 和 蛋白 连接 酶 E3 XIAP 连锁 凋 亡 抑制 蛋白 
(上 E3 ubiquitin-protein ligase XIAP，XIAP ) 基因 ,ring 
结构 域 蛋 白 1(ring-box protein 1，RBX1 ) 基因 共同 
参与 了 泛 素 介 导 的 蛋白 质 水 解 这 一 通路 ,RBXI 基 
因 表 达 上 调 , 上 调 293. 29 倍 ;XIAP 基因 表达 下 调 ， 
下 调 42.81 倍 。 




































































































































































































































































3 讨论 





人 工 诱 导 企 裙 追 寄 晶 晴 进 入 滞 育 ,可 以 有 效 延 
长 其 货架 期 。 本 研究 以 伞 裙 追 寄 蝇 非 滞 育 和 沛 育 师 
转录 组 为 研究 对 象 ,通过 对 比分 析 ,筛选 出 454 个 滞 
育 关联 基因 。 这 些 基因 主要 集中 在 信号 转 导 、 内 分 
泌 系统 、 碳 水 化 合 物 代谢 等 途径 。 

在 伞 裙 追 寄 蝇 沾 育 期 间 ,丙酮 酸 激酶 基因 丙酮 
酸 脱 氧 酶 的 El 组 分 a 亚 基 表 达 上 调 , 有 利于 中 间 
产物 一 一 乙酰 辅酶 A 生成 ,然后 进入 TCA 循环 
( Naito et al.，1998 ) 。 细 胞 色素 氧化 酶 相关 基因 的 
表达 上 调 ,利于 呼吸 电子 链 的 传递 ,为 昆虫 的 滞 育 期 
提供 能 量 。 由 此 可 以 推测 糖 酵 解 、 柠 檬 酸 循环 、 氧 化 
磷酸 化 在 企 裙 追 寄 蝇 清 育 期 间 的 能 量 供给 并 不 是 维 
持 在 一 个 低 水 平 。 在 脂 代 谢 中 ,FAR 催化 脂 酰 辅酶 
A 还 原 成 脂肪 醇 ,这 是 生成 脂肪 醇 的 主要 途径 。 关 
于 FAR 的 研究 主要 集中 在 鳞 怒 目 昆虫 信息 素 , 在 家 
百 Bombyx mori 和 梨 蛾 属 Yponomeuta 的 研究 中 表明 
一 种 FAR 催化 性 腺 中 多 种 性 信息 素 , 而 欧洲 玉米 量 
Ostrinia nubilalis 两 个 亚 种 的 研究 则 表明 2 个 FAR 
基因 生成 不 同 的 性 信息 素 (Moto et al.，2003; 
Wahlen et al.，2009 ; Lassance et al., 2010; 杨 丙 等 ， 
2012)。JHE 由 脂肪 体 细 胞 和 上 皮 细 胞 产生 ,是 主要 
的 保 幼 激素 (juvenile hormone, 理 ) 特 异性 降解 酶 之 
一 ,裂解 且 酯 键 产生 甲醇 和 JH 酸 ,在 调节 代谢 
中 发 挥 重 要 作用 ( 李 胜 等 , 2004) 。JHE 基因 表达 下 
调 , 由 此 可 以 推测 滞 育 期 间 伞 裙 追 寄 晶体 内 保 幼 激 
素 降 解 受到 阻碍 ,JH 含量 可 能 不 是 保持 在 一 个 低 
水 平 。 

肌 腿 蛋白 属于 大 细胞 外 基质 糖 蛋 白 家 族 ,在 神 
经 系统 发 育 过 程 中 以 特定 模式 表达 ,并 在 成 熟 后 下 
调 。 大 量 的 体外 研究 表明 , 肌 腿 蛋白 在 神经 系统 中 
的 作用 与 前 体 细胞 迁移 , 轴 突 生长 和 引导 有 关 
(Joester and Faissner，2001 ) 。 可 以 推测 ,在 滞 育 期 
间 例 裙 追 寄 蝇 肾 的 肌 腿 蛋白 可 能 处 于 发 育成 熟 
阶段 。 

过 氧化 氢 酶 基因 在 过 氧化 氢 物 酶 体 通路 发 挥 着 
保护 酶 的 作用 ,使 企 裙 追 寄 蝇 在 清 育 期 间 抵 抗 逆境 
胁迫 ,避免 受到 损伤 (于 德 玲 和 王 昌 留 , 2016)。 在 
菜 蛾 盘 绒 草 蜂 Cotesia vestalis 清 育 期 间 ,发 现 低 浓 度 
的 过 氧化 氢 和 高 活性 的 CAT 有 利于 增强 该 蜂 的 抗 
逆 能 力 ( 郝 仲 萍 等 , 2012; 胡 甘 雨 等 , 2014) 。 

细胞 凋 亡 是 一 种 细胞 程序 性 死亡 过 程 , 由 许多 
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过 程 相互 关联 ,共同 调控 。 在 伞 裙 追 寄 晶 沛 育 期 间 ， 
XIAP 基因 、RBXI1 基因 共同 参与 泛 素 介 导 的 和 蛋白质 
水 解 通路 , RBX1 是 泛 素 连接 酶 E3 的 一 个 催化 亚 
基 , 可 以 控制 许多 短 寿 命 调节 蛋白 的 更 新 ( Jia et al.， 
2011; 王 步 勇 等 ，2013 )。 在 秀丽 隐 杆 线虫 
Caenorhabditis elegans 的 肥 胎 发 育 ,生殖 细胞 丰 度 以 
及 成 年 活力 或 寿命 具有 重要 的 调控 作用 (Moore and 
Boyd, 2004) 。XIAP 是 具有 抑制 细胞 凋 亡 的 蛋白 ， 
RBX1 与 XIAP 发 挥 了 后 抗 作用 。RBX1 具有 泛 素 连 
接 酶 E3 活性 ,可 以 引起 XIAP 自身 泛 素 化 而 被 降 
解 ,同时 也 能 够 引起 底 物 发 生 泛 素 化 ,抑制 细胞 凋 亡 
的 线粒体 途径 ,发挥 抗 凋 亡 的 作用 ( 李 中 海 等 ， 
2003 ) 。 

本 研究 通过 对 伞 裙 追 寄 晶 非 滞 育 师 与 滞 育 晴 的 
转录 组 差异 表达 基因 KEGG 富 集 通路 的 分 析 比 较 ， 
揭示 了 其 滞 育 的 关键 基因 ,为 阐明 其 滞 育 机 理 黄 定 
了 一 定 的 理论 基础 。 关 键 滞 育 关联 基因 在 伴 裙 追 寄 
晶 清 育 过 程 中 的 差异 性 表达 及 作用 ,还 有 待 于 进 一 
步 验证 与 探究 。 
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